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Abstract:  A Streamline  Curvature  Method  of  through  flow  analysis  for 

axisymmetric  axial,  mixed  and  radial  flow . turbomachinery 
has  been  implemented.  The  method  can  lie  used  for  either 
the  direct  or  indirect  problem  and  can  solve  open  flow 
problems  such  as  flow  through  a propeller  on  a body  of 
revolution  as  well  as  the  more  familiar  duct  flow  problems. 
Computational  refinements  allow  modeling  of  blade  row* 
chordwise  and  spanwise  loading  and  blockage  distributions 
and  selection  of  streamline  locations.  Several  examples 
of  the  abilities  of  the  analysis  are  presented. 
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angle  between  a streamline  and  a reference  line 
dummy  variable  of  integration  in  n direction 
ordinate  parallel  to  reference  line 
streamline  curvature  = (radius  of  curvature  Rjc) 
ordinate  normal  to  a streamline 
static  pressure 

angle  between  streamline  and  axis  of  rotation 

2 2 1/2 

meridional  velocity  = (u  + v ) 
fluid  density 

ordinate  along  a radial  line 
ordinate  parallel  to  a streamline 
rotor  rotational  velocity 
velocity  parallel  to  axis  of  rotation 
velocity  normal  to  axis  of  rotation 
radius 

tangential  velocity 

total  velocity  at  reference  conditions 
ordinate  parallel  to  axis  of  rotation 
ordinate  normal  to  axis  of  rotation 
coordinate  location  on  n axis 


Subscripts 

x denotes  partial  differentiation  W.R.T.  X 

y denotes  partial  differentiation  W.R.T.  Y 

i denotes  inner  boundary 

o denotes  outei  boundary 
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D denotes  a function  of  n 

£ denotes  a property  of  a particular  streamline  corresponding  to 

a value  of  £ on  r) 

1 blade  row  inlet 

2 blade  row  exit 
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Introcluc  1 1 on 

The  Streamline  Curvature  Method  (SOS)  has  been  developed  over  many 
years  to  solve  various  classes  of  axisymmetric  and  quasi-three  dimensional 
turbomachine,  through  flow  problems.  The  method  relies  on  the  ability  to 
define  accurately  the  streamlines  on  a meridional  plane  and  to  determine 
radial  and  convective  accelerations  based  on  their  geometry.  In  the  past, 
the  SCM  has  been  successfully  applied  to  axial,  flow  pumps  and  compressor/ 
turbines,  this  being  performed  with  the  equations  written  in  orthogonal  or 
intrinsic  coordinate  systems.  Except  for  the  cumbersome  quasi-three  dimensional 
analysis  (the  direct  problem) , approximations  are  often  made  which  limit  the 
accuracy  of  the  analysis,  particularly  with  regard  to  the  blade  rows.  Blade 
rows  are  usually  treated  as  thick  actuator  disks  with  chordwise  loading, 
blockage,  and  radial  body  forces  ignored. 

Computationally,  many  of  these  programs  are  inefficient  and  limited 
in  the  types  of  cases  they  can  handle.  This  paper  documents  a computer 
analysis  which  addresses  the  problems  described.  The  program  is  capable 
of  modeling  axial,  mixed  and  radial  flows  with  multiple  blade  rows  and 
can  solve  an  approximately  direct  as  well  as  the  indirect  problem.  Blade 
to  blade  effects  are  incorporated  as  a circumferentially  averaged  radial 
body  force.  Particular  attention  to  stability  of  convergence  makes  this 
analysis  suitable  for  problems  which  are  usually  very  difficult  to  solve. 

In  order  to  improve  the  usefulness  of  the  SCM,  equations  were  written 
such  that  reference  stations  may  take  an  arbitrary  path  through  the  machine. 

This  allows  the  complete  modeling  of  the  leading  and  trailing  edges  of  the 
blade  rows.  The  option  of  intrablade  computing  stations  models  the  blade 
loading  and  blockage  distribution,  and  radial  body  forces.  The  use  of 
specialized  curve  fitting  routines  allows  radial  as  well  as  mixed  flow 
and  axial  flows  to  be  analyzed.  The  stability  of  the  program  allows  very 
high  station  aspect  ratios  to  be  used.  Pseudo-streamlines  and  the  ability 
to  define  streamline  locations  improves  the  accuracy  and  speed  of  the  program 
and  allows  concentration  of  data  in  regions  of  interest. 

To  date,  the  program  has  been  used  to  analyze  axial  flow  and  radial 
flow  pumps  and  has  performed  the  direct  analysis  of  an  open  propeller  in 
an  infinite  medium,  all  with  good  success.  Examples  of  these  cases  will 
be  presented  in  this  report.  Documentation  of  the  actual  program  will 
appear  in  a later  report. 

Development  of  the  Equat ions  o f Motion 

In  the  axisymmetric  inviscid  analysis,  Euler's  momentum  equation  is 
recast  into  a radial  equilibrium  equation.  Conservation  of  mass,  total 
energy  and  angular  momentum  complete  the  analysis.  The  equations  must 
be  solved  by  successive  approximations,  with  necessary  data  and  derivatives 
taken  from  the  iteratively  approximated  velocity  field  and  streamline 
geometry. 

Let  x and  y be  rectangular  coordinate's  in  a iiuiidion.il  plane  as 
illustrated  in  Figure  1.  Euler's  equations  in  these  rectangular 
coordinates,  for  a tvo-dimt  u.sional  incomp?  i . ;s ible  J low  a>  e: 


w 
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Figures  (1)  and  (2)  represent  the  meridional  plane  of  the  solution’ 
and  geometric  quantities  used  in  the  following  equations  are  shown  on 
these  figures. 


by: 


The  pressure  gradient  between  points  a and  b.  Figure  2,  is  represented 


dJP 

dn 


3F  , 

= X-  ds 
ds 


. DP  , 

+ x—  dn 
3n 


(2) 


To 


tives 


the  pressure  gradient  along  n requires  the  partial  deriva- 
determined . 


Lrcamwise  direction  one  has 


DP 

Ds 


DP 

9y 


sin 


DP 

Dx 


cos 


4> 


(3) 


Combining  (1)  and  (3)  we  have 


1 DP 
p Ds 


Note  that 


(4) 


v = V 

sin 

<J> 

and 

in 

(5) 

u = Vm 

cos 

4- 

• 

The  following  derivatives  are  determined  by  the  chain  rule: 
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_3v 

8x 


V cos  cb  d>  + sin  d>  V„, 
m ^ x mx 


8u 

3x 


-Vm  sin  <j>  <Px  + cos  <j>  vmx 


3v 

3y 


V cos  <b  (t  + sin  d>  V„, 

m Y Ty  1 niy 


3u 

3y 


-V 


m 


sin  d)  <b 

Yy 


+ cos  4)  Vmy 


(6) 


The  quantities  determined  in  (5)  and  (6)  are  substituted  into 
Equation  (4)  and  after  reduction  the  result  is: 


1 3P 
p 3s 


V cos  <j)  + V 

in  T mx  m 


sin  $ Vniy 


(7) 


Equation  (7)  reduces  to 


1 3P 

p 8s 


V 

m 


(8) 


which  is  seen  to  be  the  differential  form  of  the  steady  Bernoulli's 
equation  for  an  incompressible  flow. 

The  streamwise  normal  component  of  Equation  (2)  is  developed  in  a 
similar  manner  by  noting  that 


3P 

3n 


3P 

8y 


COS  <j) 


3P 

3x 


sin  (J> 


and 


X 3P  = 
p 3n 


3v 

3x 


+ v 


3v, 

3yJ 


3u 


3 un 


. I V U . t'  -1  . . 

COS  <J>  - [u  J-  + v ] Sin  <t> 


3y 


(9) 


(10) 


After  combining  (5),  (6)  and  (10),  we  obtain  the  following  result: 
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1 dP 
p dn 


2 2 

V cos  Ad)  + V sin  <J>  <p 
m vx  m y 


(ID 


which  is  equivalent  to 


1 

p Sn  I ds  x ds  ' y 


The  quantity  in  the.  brackets  in  (12)  is  recognized  to  be  equivalent  to 
8<{>/ds,  the  curvature  (k)  of  the  streamline. 

Finally,  Equation  (12)  becomes 


A?-  = v 2 


iDx 


+ U d» 


(12) 


i|£=kv2 

p dn  m 


(13) 


Combining  Equations  (2),  (8)  and  (13),  we  determine  the  radial  equilibrium 
equation  for  nonswiiling  flow  to  be 

dP  _ 3j?  d_s  8_P  dn 

dn  ds  dll  dn  dq 


and  finally, 


1 dP 
P 3n 


2 3Vm 

k V sin  a + V — cos  a 
m mas 


(1  A) 


The  pressure  gradient  due  to  swirl  is  determined  in  a similar  manner  and 
leads  to  the  term, 


I IP.  „ A v 2 

p dR  R 0 


(15) 
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whore  there  are.  no  circumferential  derivatives,  i.e. , the  flow  is  assumed  to 
be  uniform  in  the  circumferential  direction.  The  computational  form  of  the 
radial  equilibrium  equation  is 


I SP 
p Sn 


1 

sin  a + ~ 


9V 

sin(4>+a)  + V — cos  a 
m 9s 


(16) 


The  three  terms  on  the  right  hand  side  of  Equation  (16)  ■ epresent 
the  meridional  curvature,  the  radial  and  the  convective  accelerations 
respectively.  Integration  of  the  equation  will  yield  the  static  pressure 
difference  between  any  two  points  in  the  flow  field. 


Application  of  t he  Continuity  and  Ener gy  Equa t ion 

From  Equation  (16)  , the  static  pressure  difference  relative  to  oOise 
point,  say  n_.  , in  the  flow,  to  another  point,  along  the  reference  Dine 
may  be  found.  To  satisfy  conservation  of  mass  and  energy  across  tin: 
reference  line,  an  absolute  value  of  static  pressure  must  be  found  at  the 
reference  point,  r^.  The  following  direct  solution  for  this  value,  is  an 
improvement  over  other  procedures  which  are  iterative  in  nature. 

A continuity  equation  may  be  written  for  every  station  in  the  problem: 


V (n) 

m 


r(n) 


sin(ar|) 


dr)  = const 


(17) 


In  words,  the  mass  flow  across  a reference  line  regardless  of  its  path 
between  two  streamlines,  is  constant,  assuming  an  incompressible  fluid. 
The  energy  equation  may  be  written  for  a particular  point,  C,  on  the 
reference  line  (along  a streamline)  as 


] „ 2 

1 „ 2 

_ I 

9P 

p 4-  — o V 
«°r  2 P “V 

•* 

" 2 p 8C 

' ^ ^ ~ . 
n 

p 

i 

9>i, 

dn 

i 

P V. 


. (18) 


The  first  two  terms  on  the  left  hand  side  represent  the  total  energy 

(static  plus  dynamic)  available  at  the  point,  f, . The  term  Pn,  contains  all  head 

loss  and  rotor  energy  changes  between  the  reference  condition  and  the 

station  of  interest.  The  third  term  is  derived  from  the  conservation  ot  angular 

momentum  and  it  is  the  rotational  kinetic  energy.  The  quantity  KVg  is 

assumed  constant  in  the  absence  of  losses  or  momentum  changes  due  to  a 

rotor  or  stator.  The  sum  of  the  fourth  and  fifth  terms  is  the  static 

pressure  at  £•  The  term  in  the  integral  is  the  static  pressure,  difference 

from  Equation  (16). 
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Equations  (17)  and  (18)  may  be  combined  and  integrated: 


1 „ 2 1 „ 2 


\+Jp  \ ~ I p v°ri  - \ - 1“  p air  rn  sin(an>  dn  = 


2 P vmn  rn  sin(an)  dn 


because  is  a constant,  Equation  (19)  may  be  rearranged  to  give 


, ,1  „ 2 1 „ 2 
“n  2p\  -2pV0n 


PHi  W dn  ' 2 P V9n  - -olR- 


1 „ 2 1 tt2  , 

2 p Vmn  I v'  dn 


and  finally; 


W = rn  sin  (an) 


„ , 1 2 „ 2.  1 „ 2 

Si  + 2 p (V“r)  Rn  } - I p V, 


- dpi  w2d 
p 33  n 


The  static  pressure  anywaere  along  a reference  line  is  then: 
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pr  = pHi  + 


1 3T 

[ p 9nJ 


dn 


(22) 


T1i 


The  velocity  profile  which  satisfies  angular  momentum,  continuity  and 
total  energy  is  given  hy  rearranging  (18) : 


P 


p + o P(v 


- v6c  ) - pn.  _ 


K 

_ I 3P.' 

I 

r p 

dn 


ni 


1/2 


(23) 


The  flow  field  is  solved  by  marching  downstream  to  each  station  in  turn  and 
integrating  Equation  (16) . The  static  pressure  is  then  obtained  from  Equation 
(22).  The  improved  velocity  profiles  are  generated  by  Equation  (23)  until 
the  changes  in  the  profiles  are  small  between  tv.’o  successive  passes. 


Computati o rial  Procedure 

Certain  data  are  necessary  to  start  the  iterative  computation  cycle. 
Data  specifying  the  design,  i.e.,  the  geometry,  blade  location,  loading 
and  thickness,  and  reference  fluid  conditions  are  input.  Initial  one- 
dimensional approximations  of  all  the  velocity  profiles  and  the  streamline 
pattern  are  made.  From  the  initial  guess  for  the  velocity  profiles  and 
streamline  geometry,  derivatives  necessary  for  the  terms  in  Equation  (16) 
can  be  determined.  This  equation  is  integrated  and  the  pressure  difference 
as  a function  of  n is  saved.  This  information  is  used  to  solve  Equation 
(21)  for  Pn  . Once  solved,  Equation  (23)  is  used  to  generate  improved 
velocity  profiles.  The  improved  velocity  profiles  are  integrated  to  give 
the  mass  flow  as  a function  of  n and  by  specifying  percentages  of  the  total 
mass  flow,  new  streamline  locations  are  determined.  This  data  is  fed  back 
into  the  program  and  the  cycle  is  repeated  until  convergence  is  met. 


All  of  the  major  program  variables  such  as  streamline  location,  radius 
of  curvature  at  every  point,  velocities  and  other  derivatives  are  sever] y 
damped  against  their  previous  values  to  prevent  instabilities  from  growing 
during  the  computations.  As  a result  of  this  treatment,  problems  with  the 
ratio  of  radial  distance  to  axial  spacing  (station  aspect  ratio)  that  are 
large  can  be  solved.  This  feature  is  important  when  problems  sucli  as 
propellers  in  an  approximately  infinite  medium  are  to  be  analyzed.  The 
computational  aspects  of  the  SCM  analysis  will  he  fully  reported  in  a later 
document,  along  with  a complete  problem  solution. 
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A rotor  changes  tlic  total  head  in  proportion  to  the  change  in  angular 
momentum  of  the  fluid  passing  through  the  rotor.  The  total  pressure  is 
increased  (or  in  the  case  of  a turbine,  decreased)  in  accordance  with  the 
equation , 

AP~n  = 2 P(U2  V02  " ui  VGi>  » (24) 


where  Vq  is  positive  in  the  direction  of  rotation. 

Losses 

Frictional  and  secondary  flow  losses  play  an  important  role  in  deter- 
mining the  performance  of  any  fluid  handling  machine.  Our  ability  to  predict 
these  losses  theoretically  for  a generalized  turbomachine  in  nonexistent, 
therefore  correlational  data  must  be  used  if  the  effects  are  to  be  included 
in  the  analysis.  To  obtain  these  data,  a machine  must  be  analyzed  for 
inviscid  flow.  Comparison  with  experimental  results  then  indicates  the 
magnitude  and  distribution  of  the  losses  encountered.  These  data  can  be 
incorportated  into  the  analysis  as  a distributed  total  pressure  loss  and 
the  program  'tuned'  for  a particular  machine.  It  must  then  be  assumed 
that  the  losses  are  similar  for  other  machines  of  the  same  general  type. 

Data  exist  for  a variety  of  pumpjet  configurations  and  have  been  used 
successfully  in  predicting  the  performance  of  new  machines. 

Examples 

In  this  section  we  present  several  examples  of  the  uses  of  the  Streamline 
Curvature  Method  as  described  in  this  report.  The  first  example  is  of  a 
counterrotating  set  of  open  propellers  on  a body  of  revolution.  The  outer 
boundary  streamline  is  defined  by  a potential  flow  around  the  body  shape. 

The  effect  of  the  rotors  on  this  streamline  are  small  enough  to  be  neglected. 

The  plotted  streamlines  for  this  configuration  are  shown  in  Figures  (3)  and 
(4)  demonstrate  the  more  significant  results.  First,  there  is  a streamline 
contraction  through  the  rotors  as  the  flow  is  accelerated  near  the  body. 

The  velocity  profiles  exhibit  the  characteristic  bulge  or  jet  behind  the 
rotor  and  in  this  case  the  average  jet  velocity  is  about  1.6  times  the  free 
stream  velocity.  The  plot  also  demonstrates  the  ability  of  the  program 
to  use  curved  reference  stations  and  to  solve  problems  with  high  station 
aspect  ratios  (in  this  case,  AR=25) . The  analysis  is  able  to  determine 
either  the  tip  radius  of  the  propellers  for  a given  mass  flow  rate,  or, 
given  the  tip  radius  as  in  the  direct  solution,  to  determine  the  mass  flow 
and  powering  requirements  for  the  configuration. 

The  second  example  is  of  a Francis-type  turbine  with  wicket  gates  and 
a simulated  inlet  volute.  This  particular  problem  demonstrates  the  ability 
of  the  analysis  to  describe  axial,  uvi\<d  and  purely  radial  flows.  A 
streamline  plot  is  presented  in  Figure  5. 
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The  third  example  is  the  direct  solution  of  the  flow  through  an  open 
propeller . In  this  case  only  the  angularity  distribution  in  the  rotor 
exit  plane  is  specified,  rather  than  the  usual  tangential  velocity  distribution. 
A tangential  velocity  distribution  is  iteratively  determined  that  satisfies 
the  angularity  while  at  the  same  time  all  other  equations  of  motion  are 
satisfied.  Figure  (6)  shows  t lie  theoretical  velocity  profile  obtained  from 
this  procedure  and  corresponding  experimental  data  for  one  typical  case. 


Summary  ' - ' 

A method  of  through-flow  analysis  for  turbomachinery  has  been  developed 
which  lias  proven  successful  for  a variety  of  types,  including  axial,  radial 
and  mixed  flow  machines.  The  equations  of  motion  are  solved  by  a different 
approach  than  is  commonly  used,  in  that  a differential  equation  for  the 
pressure  gradient  rather  than  the  velocity  gradient  is  used.  Additionally, 
the  energy  and  continuity  equations  are  satisfied  by  a direct  integration 
rather  than  by  iterative  approximations. 

Refinements  to  the  numerical  method  permit  more  accurate  modeling  of 
the  streamlines,  allowing  mixed  and  radial  flows  to  be  analyzed.  Multiple- 
blade  rows  with  intrablade  computing  stations  allows  modeling  of  the  blade 
spanwise  and  chordwise  loading  distributions  and  the  blade  thickness 
distribution.  Selection  of  streamline  location  and  density  improves  accuracy 
in  regions  of  interest. 

Particular  attention  to  the  stability  of  the  numerical  procedure  allow:; 
very  high  station  aspect  ratios  to  be  used  and  allows  problems  such  as 
propellers  in  an  unbounded  medium  to  be  analyzed,  as  well  as  the  standard 
internal  flow  problems. 

Friction  losses  and  secondary  flow  effects  may  be  included  in  an 
empirical  fashion  when  a particular  machine  type  is  prescribed. 


r 
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Figure  3 - Overall  View  of  the  Streamlines  Through 
a Counterrotating  Open  Propeller  Set 


-19-  21  July  1977 

MWM: jep 


DISTANCE  FROM  PROPELLER  HUB  (inches) 


Figure  6 


-21- 


21  July  1977 
MWM: jep 
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Comparison  of  Theoretical  and  Experimental 
Velocities  Behind  an  Open  Propeller  as  Solved 
by  the  Direct  (Specified  Angularity)  Method 
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